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Abstract 

Previous studies of the QED systematic uncertainties on the LEP measurement of the W-boson 
mass, Mw, have used idealized event selections and fitting procedures. In this paper, the Monte 
Carlo tandem of KoralW and YFSWW is used to investigate how the full experimental Mw 
extraction procedure affects these estimates. It is found that the kinematic fitting used in 
the experimental determination of Mw enhances the sensitivity to QED corrections involving 
real photon production. It is concluded that the previous estimates of the QED theoretical 
uncertainty on the LEP2 Mw measurement may be too small. A simple procedure for approx- 
imating the effect of the kinematic fit at the level of the generated four fermions is proposed. 
This procedure would allow previous theoretical studies to be repeated using a much closer 
approximation of the real experimental mass extraction method. Finally, the possibility of 
setting experimental limits on 0(a) theoretical uncertainties using identified e + e~ — ► W + W~7 
events is discussed. 
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1 Introduction 



It is anticipated that the final experimental error on the mass of the W boson, M w , from LEP 2 
will be approximately 35 MeV. A potential source of theoretical uncertainty is the treatment 
of QED and electroweak corrections in the Monte Carlo programs used to simulate the process 
e + e~ — > 4/(7). In recent years 0(a) electroweak corrections have been implemented in the 
YFSWWfl] and RacoonWW[2] programs. The current estimate of the theoretical systematic 
uncertainty on the LEP measurement of M w due to these predominantly QED corrections is less 
than 5 MeV [3]. However, as pointed out by the authors, this estimate is based on an idealized 
M w extraction procedure. Indeed, the authors suggest that the study should be repeated using 
the full experimental fitting procedure. 

In this paper the concurrent Monte Carlo combination [4] of KoralW[5] and YFSWW[1] 
(referred to as KandY) is used to study the sensitivity of the full experimental analysis to 0(a) 
corrections. Monte Carlo events generated using KandY are passed through the OPAL detector 
simulation[6] and W-boson reconstruction algorithm[7]. In this way, the previous estimates 
of the QED systematic uncertainties can be extended to the full experimental Mw extraction 
procedure. It is found that the sensitivity to the so-called[3] 0(a) non-leading (NL) electroweak 
corrections is significantly enhanced by the kinematic fitting used in the experimental analyses. 
This enhancement is a result of the relatively large change in the cross section for production 
of e + e~ — > W + W~7 events where the photon is away from the region collinear with the e + e~ 
beams. Consequently, measurements of the cross section for e + e~ — > W + W~7 in this region 
place constraints on the related M w systematic uncertainties. 

It is suggested that previous studies [3] be repeated using a closer approximation of the real 
experimental procedure. With this in mind, a simple procedure to approximate the effects of 
the kinematic fit at the level of the generated four fermions is proposed. 

2 The Experimental Measurement of -M\y at LEP 2 

The main measurements of Mw at LEP are obtained from the direct reconstruction of the 
W-boson invariant mass distribution in e + e~ — > qqiVe and e + e~ — > qqqq events. Although the 
exact methodologies used by the four LEP experiments differ [7, 8], there are many common 
features and the general principles of the analyses are summarized below. 

The experimental event selections for e + e~ — > W + W~ have high efficiency, typically 80 % — 
85 %. Once events are selected, algorithms are applied to reconstruct the four- momenta of 
the fermions. In e + e~ — > qqiVf events, the charged particle tracks and any clusters in the 
calorimeters (electromagnetic and hadronic) which are not associated with the lepton are forced 
into two jets using the Durham[9] or LUCLUS[10] algorithms. In e + e~ — > qqqq events, tracks 
and clusters are forced into four jets (although DELPHI and OPAL allow for the possibility of a 
fifth 'gluon' jet). The jet energy resolution is relatively poor, &e/E > 60 %/y/E, and neutrinos 
from leptonic W-boson decays are unobserved. Consequently, kinematic fits are employed to 
improve the event-by-event mass resolution. The kinematic fits impose the constraints of energy 
and momentum conservation. In most analyses the additional constraint of equal masses for 
the two W bosons is imposed, effectively neglecting Tw, referred to as a five constraint (5C) 
fit. The result of the 5C kinematic fit is a single reconstructed mass for each event which 
approximates to the average of the invariant masses of the fermion pairs from the two W-boson 
decays. In general, photons from initial state radiation (ISR) within the detector acceptance 
are not treated as separate particles in the fit. Although the four LEP collaborations employ 
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different procedures to extract from the reconstructed mass distribution, each method 
essentially corresponds to fitting the peak of the reconstructed mass distribution and calibrating 
out experimental biases using Monte Carlo event samples. 

The single most important aspect of the kinematic fit is the constraint that the sum of 
the energies of the four fermions is equal to the centre-of-mass energy of the e + e~ collision, 
y/s. Due to the poor jet energy resolution, this constraint significantly improves the invariant 
mass resolution. However, in the presence of ISR, this procedure introduces a bias in the 
reconstructed mass, since the energies of the four fermions should be constrained to the centre- 
of-mass energy after photon radiation, y/s', rather than to y/s. As a result, for events with ISR, 
the kinematic fit will return a mass value which is too large. On average, the event-by-event 
mass returned by the kinematic fit approximately corresponds to 



Hence, the reconstructed invariant mass distribution depends strongly on the distribution 
and as a result the peak of reconstructed W-boson mass is approximately 100 MeV higher than 
M w (the exact value depending on y/s). The bias in the fitted mass is even greater, typically a 
few hundred MeV. This bias is subsequently determined using Monte Carlo events and the data 
measurement corrected. Consequently, deficiencies in the Monte Carlo simulation may result 
in systematic biases and it is necessary for the Monte Carlo programs to predict accurately the 
xTs* distribution. 

3 QED Uncertainties on the LEP M" w Measurement 

The current estimate of 5 MeV for the QED/Electroweak theoretical uncertainty on the LEP 
Mw measurement is based on the Monte Carlo generated invariant mass distribution of the 
/i~z7 M system in e + e~ — * fj,~V ^ud^) events[3]. Uncertainties are obtained from the shifts in 
the invariant mass distribution which result from dropping higher order theoretical terms. The 
simplest definition of invariant mass is obtained directly from the four-momenta of the gen- 
erated fermions (termed BA RE mass 1 ). However, in the real experimental analyses, photons 
from final state radiation (FSR) are usually recombined with the nearest jet or lepton. This 
recombination is emulated at the level of the generated particles by recombining photons with 
the closest fermion if the photon has energy E 1 < 1 GeV or if the invariant mass of the photon 
and fermion is less than 5GeV(25GeV), termed the CAL05(CAL025) photon recombination 
scheme. In this way FSR photons are usually recombined with the appropriate fermion. In 
the CAL05/CAL025 schemes photons within 5° of the beam direction are treated as unob- 
served and are not recombined, reflecting the holes in experimental acceptance along the beam 
direction. 

The main difference between the definitions of BARE and CALO masses and the ex- 
perimentally reconstructed mass arises from the kinematic fit employed in the experimental 
analyses. As a result of the kinematic fit, 0(a) QED corrections may affect the reconstructed 
mass in two ways: either by modifying the invariant mass distribution of the fermion pairs 
from W decays or by modifying the \/iP distribution. In the first case, the theoretical estimates 
of mass biases obtained using either BARE or CALO schemes provide a good estimate of the 
true experimental bias. However, in second case, even if the invariant mass distribution of the 

1 The nomenclature of reference [3] is used throughout this paper rather than that of [11]. 
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fermion pairs is unchanged, the experimentally measured mass distribution may be biased as a 
result of the kinematic fit scaling the measured energies to y/s rather than \/P. Consequently, 
previous estimates of the impact of theoretical corrections which modify either the rate or 
energy spectrum of real photons produced in the process e + e~ — > W + W~7 may be too small. 

To study possible differences between the previous methods used to estimate systematic 
uncertainties and the true impact on the experimental measurement the Monte Carlo tandem 
of KORALW and YFSWW (KandY) is used. A total of four million e+e" -> W+W" events 
was generated over 8 separate values of y/s reflecting the LEP 2 data sample used in the W- 
mass determination (183 GeV < y/s < 209 GeV). The events were passed through the full 
OPAL detector simulation[6] and W-mass reconstruction procedure including the 5C kinematic 
fit [12]. The following studies are based on events passing the OPAL W + W~ — > qqiVe and 
W + W~ — > qqqq event selections. The KandY program produces correction weights enabling 
generated events to be reweighted to the corresponding theoretical prediction removing certain 
terms. Correction weights were produced corresponding to: the replacement of the Screened 
Coulomb correction [13] with the unscreened correction; degrading the treatment of ISR from 
0(a 3 ) to 0(a 2 ) YFS exponentiated [14] leading logarithm (LL); and dropping the so-called 
0(a) NL electroweak corrections of KandY. It is important to note that the corrections termed 
l O(a) NL' in [3] include bremsstrahlung from the W bosons (referred to as WSR) and the 
resulting interference with ISR. 

The treatment of the Coulomb correction modifies the invariant mass distribution of the 
fermion pairs from the W decays. Consequently, the shift in the BARE jCALO invariant mass 
distribution provides a good measure of the true experimental bias. Degrading the treatment 
of ISR from 0(a 3 ) to 0(a 2 ) exponentiated LL has negligible impact on the experimental 
measurement of Mw [7] (the measurement is sensitive to ISR through a modification of the x/s 1 
distribution, however the effects of the 0(a 3 ) terms are small). The effect of the 0(a) NL 
electroweak corrections, which include bremsstrahlung from the W bosons, is more interesting. 
To study potential changes in distributions it is convenient to define the fractional change in a 
differential cross section, 

(da^t - da KY ) 

^NL = j , \£) 

aa KY 

where daxy is the differential cross section in a given bin from the full KandY generation and 
da^Y 1S the corresponding distribution without the 0(a) NL electroweak corrections. 

Figure 1 shows the effect of the O(a) NL corrections on the photon energy and polar angle 
distribution. These plots are similar to those of Reference [15] showing that the experimental 
selections do not introduce significant biases for events with additional photons 2 . The inclusion 
of the O(o) NL corrections reduces the number of real photons produced in the detector ac- 
ceptance. The maximum fractional change occurs in the region perpendicular to the beam axis 
where the W + W~7 cross section is decreased by 30 %. This modification of the photon rate 
and angular distribution is predominantly due the inclusion of radiation from the W bosons. 
Specifically, the reduction in the W + W"7 cross section results from interference between ISR 
and radiation from the W bosons 3 . The photon energy spectrum is not greatly distorted by 

2 The plots in Reference [15] are for W + W~ — > qq£i>e7 events generated at y/s = 200 GeV and do not include 
the effects of event selection. 

3 This was verified using YFSWW. The reduction in the cross section for photon production shown in Fig- 
ure la was reproduced by switching from KEYC0R=2 to KEYC0R=3 in YFSWW, corresponding to switching from 
YFS exponentiation for ISR alone to the full YFS form factor including radiation from the W bosons (WSR) 
and interference between ISR and WSR. 



3 



the NL corrections as can be seen from Figure 1. 

The potential bias in the measurement of M w due to the 0(a) NL terms is investigated 
using four definitions of reconstructed mass were considered: 

• M bare '■ the average of the two W-boson masses determined from the four- momenta of 
the four fermions. 

• McALOb '■ the average of the two W-boson masses determined from the four-momenta of 
the four fermions after applying the CALOb photon recombination procedure. 

• M 5 c '■ the mass returned by the OPAL 5C kinematic fit. Only events for which the fit 
converges with a probability of greater than 0.1 % are plotted. 



• M C alo5* '■ M C alo5 rescaled by . / -4 where ^/s 1 is the invariant mass of the four fermion 



The CAL05* definition, proposed in this paper, represents an attempt to emulate the 
effects of the kinematic fit at the level of the four-momenta of the fermions. The distortions of 
the W-mass peak, plotted for above definitions of invariant mass, are shown in Figure 2. To 
quantify the effect on the W-mass measurement, the resulting biases in the mean reconstructed 
W-boson mass (in the region 75 GeV — 90 GeV) are listed in Table 1. The inclusion of the 
0(a) NL corrections has almost no effect on the BARE mass distribution and the mean value 
changes by just 1 MeV. A small but visible effect is observed for the mass in the CAL05 scheme 
giving a bias in the mean mass of 3 MeV. These results are consistent with those presented 
in [3]. However, the inclusion of the 0(a) NL corrections produces a significant shift in the 
experimentally reconstructed mass from the 5C kinematic fit, with the mean reconstructed mass 
in the range 75 GeV — 90 GeV shifted by 21 MeV. Although the invariant mass distribution of 
the fermion pairs is almost unaffected the reconstructed mass distribution is significantly biased. 
The origin of this modification of the x/s' distribution is the reduction in the production cross 
section for e + e~ — > W + W~7, which results from the interference of ISR and radiation from 
the W bosons. It should be noted that by scaling the mean mass in the CAL05 scheme by 
\/s/s' (denoted CAL05*) the effect of the kinematic fit can be approximated at the level of 
the generated fermions. 

The effect of the kinematic fit is to make the experimental analyses more sensitive to 
theoretical corrections which modify the rate/spectrum of photons produced in the process 
e + e~ — >■ W + W~7. The bias of 21 MeV associated with the 0(a) NL corrections in KandY is 
almost an order of magnitude larger than previously quoted. Since the origin of the 21 MeV 
bias is interference between ISR and radiation from the W bosons, one might expect this to be 
well accounted for in YFSWW, but this requires further evaluation. It should be noted that 
the 21 MeV bias would be reduced significantly if observed photons were treated as additional 
particles in the kinematic fits. In this case the energies of the fermions from W decay would 
be scaled correctly to v^. 

For the OPAL W-mass extraction procedure the effect of the O(a) NL corrections on the full 
Mw reconstruction procedure is somewhat smaller than the shift in the average event-by-event 
reconstructed mass quoted above. Using a modified version of the OPAL M w convolution fit [7] 
the bias associated with dropping the O(a) NL corrections is 14 MeV. 
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Scheme 


AM NL 


BARE 
CAL05 
CAL05* 
5C FIT 


+ lMeV 
+3 MeV 
+22 MeV 
+21 MeV 



Table 1: The change in the mean W-boson mass in the range 75 GeV — 90 GeV when dropping 
the electroweak 0(a) NL corrections from KandY. Shifts are shown for three different theoreti- 
cal definitions of the event-by-event reconstructed mass: BARE, CAL05 and the new definition 
proposed in this paper, CAL05*. 5C FIT refers to the shift in the average of reconstructed 
mass distribution for the full OPAL event reconstruction using a 5C kinematic fit. 

4 Limits from the Measurement of the e + e _ — > W + W~7 
Cross Section 

Measurements of the process e + e~ — > W + W~7 have been performed at LEP[16, 17, 18]. These 
measurements, which are restricted to photons within the experimental acceptance (typically 
|cos# 7 | < 0.95), provide a direct probe of 0(a) QED corrections which affect real photon 
production away from the region collinear with the e + e~ beams. Of particular interest here are 
the 0(a) NL electroweak corrections of KandY which decrease the calculated e + e~ — > W + W~7 
cross section. The largest change in the differential W + W~7 cross section occurs close to the 
beam direction (as would be expected from ISR-WSR interference). However, when integrated 
over the photon polar angle, 88 % of the change in cross section occurs within the nominal 
experimental acceptance of |cos# 7 | < 0.95. Hence, by measuring the e + e~ — > W + W~7 cross 
section for events where the photon is observed, it is possible to probe experimentally the 0(a) 
NL corrections of KandY. 

When the four LEP experiments publish final results including all LEP 2 data, the com- 
bined experimental precision on the W + W~7 cross section should be approximately 7 %. The 
inclusion of the O(o) NL electroweak corrections of KandY produces a 20 % decrease in the 
e + e~ — > W + W~7 cross section for |cos# 7 | < 0.95 and E 1 > 5 GeV (cuts chosen to match the 
experimental acceptances of the LEP detectors). Thus, using rate alone, the measurements are 
of sufficient statistical precision to provide a useful test of the KandY implementation of WSR 
and in particular of WSR— ISR interference. Agreement with the predictions of KandY would 
provide the first indirect observation (albeit only at the 3 a level) of real photon radiation from 
the W bosons through interference with ISR diagrams. The statistical precision of this test 
would be improved slightly by including angular information in addition to the photon rate. 

In the context of the measurement at LEP the inclusion of the 0(a) NL corrections 
in KandY produces a 21 MeV shift in the mean reconstructed mass which results from the 
modification of the \fs' distribution. This is a direct consequence of the decrease in the e + e" — > 
W + W~7 cross section. Therefore, measurements of the e + e~ — > W + W~7 cross section at LEP 
provide constraints on the 0(a) QED uncertainties related to real photon production away from 
the region collinear with the beams. A 20% reduction in the W + W"7 cross section corresponds 
to a 21 MeV shift in the reconstructed value of M\y. Assuming the bias scales linearly, a 
measurement of the W + W~7 cross section with a statistical precision of 7 % would correspond 
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to approximately 7MeV on the associated systematic uncertainty on M\y. It should be noted 
that this limit would apply to all theoretical corrections to the process e + e~ — > W + W~ which 
modify the rate of real photon production away from the collinear region; it is not specific to 
the 0(a) NL corrections of KandY. Given that the expected statistical uncertainty on the LEP 
M\v measurement is approximately 35MeV, a LEP combined measurement of the W + W~7 
cross section would provide a useful experimental constraint on the 0(a) QED systematic 
uncertainty. 

5 Conclusions 

Previous studies of QED systematic uncertainties on the LEP measurement of M w have used 
idealized event selections and fitting procedures. In this paper the KandY Monte Carlo pro- 
gram was used to study the effect of the experimental methodology on the estimated QED 
systematic uncertainty. It is found that the kinematic fitting used in the experimental determi- 
nation of Mw enhances the sensitivity to all QED corrections involving real photon production 
and it is concluded that previous estimates of the QED theoretical uncertainty on the LEP 2 
M w measurement may be too small. A simple procedure for approximating the effect of the 
kinematic fit at the level of the generated four fermions is proposed. This procedure would 
allow previous theoretical studies to be extended to include a good approximation of the real 
experimental mass extraction method. Finally, it is shown that useful experimental limits on 
O (a) theoretical uncertainties related to real photon production away from region collinear 
with the e + e~ beams may be inferred from LEP measurements of the e + e~ — > W + W~7 cross 
section. 
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Figure 1: The effect of the 0{a) NL corrections in KandY on the distributions of the energy 
and angle with respect to the e~ beam of the highest energy photon in the event. Only photons 
which have not been recombined using CAL05 scheme are included. 
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Figure 2: Fractional change in the differential cross-section as a function of the reconstructed 
mass using the BARE, CAL05, and CAL05* schemes compared to the experimental bias 
using the OPAL 5C kinematic fit. 
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